Abstract

A new class of bivariate distributions is deduced by specifying its conditional hazard functions (hfs) which are Kumaraswamy distribution. The interest of this model is positively, negatively, or zero correlated. Properties and local measures of dependence of the bivariate Kumaraswamy conditional hazard (BKCH) distribution are studied. The estimation of type parameters is considered by used the maximum likelihood and pseudolikelihood of the new class. A simulation study was performed to inspect the bias and mean squared error of the maximum likelihood estimators. Finally, an application is obtained to clarify our results with the maximum likelihood and pseudolikelihood. Also, the results are used to compare BKCH distribution with bivariate exponential conditionals (BEC) and bivariate Lindley conditionals hazard (BLCH) distributions.

1. Introduction

The problem of constructing new families of continuous bivariate distributions is one of the important and continue active research problem in statistics. This is due from one side to the limitations of the existing distributions and its lack in modeling some stochastic phenomena and from the other side to meet the needs of modern and developed applied field. Balakrishnon et al. [1] obtained a new technology to build the bivariate probability density functions (pdfs) with specified conditional hfs, and the development of the models determines many of the several properties, including survival function (sf), hazard bivariate function, the Clayton–Oakes measure, and the measure conditional pdfs and hfs of the marginal and conditional pdfs. Moreover, Navarro and Sarabia [2] studied the reliability properties for the series and parallel systems with component lifetimes having these dependence models.

Therefore, the authors introduce a new form of the bivariate distributions with Kumaraswamy conditional hfs, and we studied some of the characteristics and discussed the estimation procedures of the parameters of this distribution. Moreover, we estimate the parameters of this distribution by two methods, e.g., maximum likelihood estimation (MLE) and maximum pseudolikelihood estimation (MPLE) procedures. Finally, this distribution is compared with other distributions considering a real dataset collected from a secondary source (p. 374 of Johnson and Wichern [3]).

A continuous random variable is Kumaraswamy distribution if its pdf (Kumaraswamy [4] is

The corresponding sf isand the hf is

We note that

This distribution is applicable to some scientific, practical, and hydrological applications (Dey et al. [5]; Ishaq et al. [6]).

2. BKCH Model

Suppose that the conditional hfs (Kumaraswamy) are distributed aswhere and are function of and , respectively.

Then, the bivariate distribution from (6) will be called a bivariate BKCH distribution.

Theorem 1. The general bivariate distribution with conditional hfs as and iswhere is the normalizing constant such that integrates to 1.

Proof. The conditional pdfs from (6) using (5) areThen, the identity yields the following relation:where and denote the marginal pdfs.
Denoting then taking logarithms in (9) and using (10) and (11) WITH then taking logarithms in (10), using (11) and (12)which is a functional equation of the form , whose most general solution is given by Aczel [7] asFrom (11) and (12), we have the corresponding marginal pdfs given byandEquation (7) describes the complete class of BKCH distribution that has the three parameters , and .
Figure 1 shows the BKCH given by (7) for the special cases for any , , and .
The conditional hfs for the BKCH distribution areThe compatibility of (17) and (18), Balakrishnan et al. [1] secures the existence of BKCH distribution. Figure 2 shows hf of the BKCH given by (7) for the special cases for and Similarly,

3. BKCH Properties

3.1. The Conditional Distributions

The conditional pdfs of the BKCH distribution arei.e.,

When , the conditional pdfs are identical.

and given by (19) and (20) satisfy the compatibility conditions studied by Arnold et al. [8] and secure existence of BKCH distribution.

3.2. Regression Function

Now, using (19) and (20), the conditional moments are

For we haveand

We notice that the regression function for is nonlinear increasing for and decreasing for (Figure 3). Additionally, the regression function for is nonlinear increasing for and decreasing for (Figure 4).

3.3. Stochastic Ordering

Let and be two-bivariate BKCH distribution, and such that and and increasing in x, for all Also, decreases in y, for all Based on the results obtained by Shanthikumar [9], then we obtain likelihood ratio order (, hazard rate order , and stochastic order , for all .

4. Local Measures of Dependence

Holland and Wang [10, 11] introduced a local dependence function given by the following formula:where is the joint pdf of x and y.

Theorem 2. The BKCH distribution given by (8) is TP₂ (totally positive of order 2) iff and TN₂ (totally negative of order 2) iff .

Proof. From (7), we obtainHence, . If , then BKCH is , and . If , then BKCH is .

Remark 1. (1)It follows from theorem (2) that if the BKCH distribution, given by (7), is or , then X and Y are positively or negatively, respectively, correlated.(2)The conditional pdfs in (19) and (20) have the local measures of dependence:(3) iff ; furthermore, the two random variables X and Y are independent iff Balakrishnan and Lai [12].The class has two parameters only by choosing ; the joint pdf (7) can be written asFigure 5 shows the joint pdf given by (29) for the special cases for and .
The conditional hfs for the joint pdf given by (29) areand the conditional pdfs areNow, the marginal pdfs of X and Y using (29) areAnd the moment of joint pdf given by (29) is

5. Estimation of BKCH Parameters

Suppose that are observed from BKCH distribution.

5.1. MLE of BKCH Parameters and Asymptotic Confidence Intervals

From the BKCH distribution, the log-likelihood function for the sample is

The likelihood equations are

The MLE can be obtained numerically using systems (35). For the asymptotic confidence interval (CI), the normal approximation of the MLE can be used to construct asymptotic CIs for the parameters when the sample size is large enough. A two-sided (1 − α)100% CIs for are , where are the asymptotic variances of .

5.2. MPLE of BKCH Parameters

In [13, 14], the study of alternative estimation technique is not based on normalizing constant. The MPLE for joint probability mass function (pmf) or joint pdf is based on the maximization of the product conditional pmfs or conditional pdfs, respectively. It is clear from (19) and (20). The log pseudlikelihood () function is

The MPLE of BKCH parameters can be obtained by solving the equations:

6. Simulation Study

In this section, the MLE and MPLE approaches are used to estimate the parameters , , and of the BKCH distribution. The population parameters are generated using software Mathematica package. The sampling distributions are obtained for different sample sizes n = 50, 100, 200, 300, 400, and 500 from N = 500 repetitions. This study presents an assessment of the properties for both MLE and MPLE techniques in terms of bias and mean square error (MSE). A general form to generate x from one marginal and then simulate a corresponding a bivariate vector using the conditional density is . The MLEs are reported in Table 1 for BKCH (1.5, 3.5, 0.05) and Table 2BKCH (1.5, 2.5, -0.05) with .

7. Application of Real Data

Suppose that X and Y are Dominant Ulna and Ulna bones. Then, the bivariate data in Table 3 represent the bone mineral density (BMD); after one year of birth, measure in gm/cm2 for 24 kids [3].

The statistic measures for the given data are The results in Table 4 represent the MLE and MPLE of BKCH parameters.

The joint pdf of bivariate exponential conditionals distribution is [15]where C is the normalizing constant.

The joint pdf of bivariate Lindley conditionals hazard distribution is [16]where is the normalized constant. is ; also, the conditional pdfs and are Lindley distribution with parameter and , respectively. Table 5 includes log-likelihood and AIC and BIC of BKCH with other models.

We note that AIC and BIC of the BKCH model more than the corresponding of the BLCH and BEC models which means that BKCH distribution is better to fit for the given data. The approximate 95% two-sided CI of the parameters , and are given, respectively, as [], [], and [, ].

8. Conclusion

In this study, we introduced a BKCH distribution based on specified conditional hfs of Kumaraswamy distribution and its local measures of dependence. In addition, the methods of MLE and MPLE of BKHC parameters are present. In view of results, Theorem 2, the interest of BKHC is positively, negatively, or zero correlated; this indicates the generality of the distribution. Furthermore, the simulation results showed that MLE operates quite uniformly, and it can be used to estimate the BKHC parameters. In particular, in Table 4, the MPLE is better than MLE because the MPLE technique uses the conditional pdfs which does not contain the normalizing constant. Therefore, Table 5 shows the BKCH distribution is a better fit for the given data compared to the BLCH and BEC distributions.

Data Availability

All data are available within the article.

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Acknowledgments

The authors are thankful to Taif University and Taif University researchers supporting project no. TURSP-2020/160, Taif University, Taif, Saudi Arabia.